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ABSTRACT 

We study the dependence of galaxy properties on the clustercentric radius and the environment 
attributed to the nearest neighbor galaxy using the Sloan Digital Sky Survey (SDSS) galaxies as- 
sociated with the Abell galaxy clusters. We find that there exists a characteristic scale where the 
properties of galaxies suddenly start to depend on the clustercentric radius at fixed neighbor environ- 
ment. The characteristic scale is I ~ 3 times the cluster virial radius depending on galaxy luminosity. 
Existence of the characteristic scale means that the local galaxy number density is not directly re- 
sponsible for the morphology-density relation in clusters because the local density varies smoothly 
with the clustercentric radius and has no discontinuity in general. What is really working in clusters 
is the morphology-clustercentric radius-neighbor environment relation, where the neighbor environ- 
ment means both neighbor morphology and the local mass density attributed to the neighbor. The 
morphology-density relation appears working only because of the statistical correlation between the 
nearest neighbor distance and the local galaxy number density. We find strong evidence that the hy- 
drodynamic interactions with nearby early-type galaxies is the main drive to quenching star formation 
activity of late- type galaxies in clusters. The hot cluster gas seems to play at most a minor role down 
to one tenth of the cluster virial radius. We also find that the viable mechanisms which can account 
for the clustercentric radius dependence of the structural and internal kinematics parameters are ha- 
rassment and interaction of galaxies with the cluster potential. The morphology transformation of 
the late-type galaxies in clusters seems to have taken place through both galaxy-galaxy hydrodynamic 
interactions and galaxy-cluster/galaxy-galaxy gravitational interactions. 

Subject headings: galaxies: clusters: general - galaxies: evolution - galaxies: formation - galaxies: 
general 



1. INTRODUCTION 

In galaxy clusters the average morphology of galax- 
ies changes with clustercentric radius or local density. 
The central region of clusters at the present epoch is 
dominated by early-type galaxies. This is known as 
the morphology-radius or morphology-density relation 
(hereafter MRR and MDR, respectively). Since the 
local density is on average a monotonically decreasing 
function of clustercentric ra dius, they appear to con- 
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Smit h et al . 2005; lPostman" et al . 200a iWemmam7 et alJ 



2006). However, there is discussion about which one 
plays a critical role in determining the galaxy mor- 
phology between the clustercen t ric radius and the local 
density (IWhitmore etafl 119931: iDommguez et al.l 120011 : 
iGoto et al.ll2003tlThomas fc Katgertll2006f) . There is also 
a remaining issue what physical parameters of galax- 
ies (morphology, color, star formation rate, or stel- 
lar m ass) correlate fundamentally with the environment 
(e.g.JChristlein fc Zabludofl1l2005t iQuintero et aLl l2006l; 
iPoggianti et al.ll2008HSkibba et al M2008h . 

A number of physical mechanisms have been pro- 
posed to explain this morphology-environment relation. 
It is suggested that tidal interactions between individ- 
ual galaxies can change mass profile and transfor m disk 
galaxies into spheroidals (jBvrd fc Valtonenl Il990) . Fre- 
quency of galaxy-galaxy interaction depends strongly on 
the clustercentric radius, and thus galaxy interactions 
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can result in the MRR or MDR. However, the duration of 
tidal interactions and thus the total tidal energy deposit 
are expected to be too small to change galaxy properties 
significantly for cluster galaxies that are moving at high 
speeds. Due to the fast orbital motions the frequency 
of galaxy mergers is also small in clusters. A series 
of frequent high-speed tidal interactions among cluster 
galaxies, called harassment, can produce im pulsive heat- 
ing a nd cause morphology transformation (|Moore et al.l 
1996). Tidal interaction with the whole cluster potential 
well can also change the structure and activity of g alaxies 
significantly (jMoss fc Whittle! 120001 : lGnedinll2003l ). 

Existence of the hot X-ray emitted by intracluster gas 
makes the mechanisms relying on hydrodynamic pro- 
cesses appear attractive. Stripping of cold gas from 
infalling late-type galaxies by a ram pressure of the 
hot intracluster gas, was proposed to explain the in- 
crea se of the early-type fraction toward the cluster cen- 
ter (jGunn fc Gottl[l972l ). Hot gas removal from infalling 
galaxies through hydrodynamic interactions with the in- 
tracluster medium can shut off gas supply and star for- 
mation activity (SFA) after the already existing cold gas 
is consumed. The interstellar medium of a galaxy travel- 
ing in the hot intracluster medium can be stripp ed from 
the d isk due to a viscosity momentum transfer (|Nulsenl 
1982) or evapo rate as the temperature o f the interstellar 
medium rises (jCowie fc Songailal Il977| ). These hydro- 
dynamic processes have a common drawback that they 
can turn spirals only to SO's and do not produce ellipti- 
cals. The major effect of the hydrodynamic processes is 
quenching the SFA by removing or ionizing the cold gas 
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in the disk of late-type galaxies. The bulge-to-disk ratio 
is not expected to increase by these processes. Therefore, 
none of the mechanisms proposed so far is able to fully 

account for the MRR o r MDR obse rved in cluste rs. 

Recently, iPark et all (|2008f ) and iPark fc Choi (|2009f ) 
have found that galaxy properties like morphology and 
luminosity depend strongly on the distance and morphol- 
ogy of the nearest neighbor galaxy. This dependence 
was found even when the large-scale background density 
is fixed, and thus is completely different from the com- 
monly known morphology-local density relation. Most 
importantly, the effects of the nearest neighbor change 
at the characteristic scale given by the virial radius of 
the neighbor galaxy. When a galaxy is located within 
the virial radius of its nearest neighbor, its morphology 
tends to be the same as that of the neighbor. But such 
tendency disappears when the separation is larger than 
the virial radius. This fact strongly suggests that galaxies 
interact hydrodynamically when the separation to their 
nearest neighbor is smaller than the virial radius of the 
neighbor, and that the effects of such interaction are sig- 
nificant enough to change their morphology and SFA. 

Outside the virial radius galaxy morphology still de- 
pends on the distance to the nearest neighbor, but is 
suddenly independent of neighbor's morphology. The 
probability for a galaxy to be an early type monotoni- 
cally decreases as the separation increases. This strongly 
supports the idea that the conformity in morphology 
of galaxy pairs is not primordial but an acquired one 
through interactions since there should be no reason for 
the break in morphology conformity to be at the cur- 
rent virial r adius of the neighbor if it is initially given. 
IPark et all (|2008l ) proposed that the mechanism respon- 
sible for this be tidal interactions causing both gravi- 
tational and hydrodynamic effects. Outside the virial 
radius the interactions are purely gravitational between 
the galaxy plus dark halo systems, but within the virial 
radius hydrodynamic (and radiative) effects must be in- 
volved as well. 

A scries of tidal interactions and mergers will keep 
the total cold gas contents in bright galaxies decreas- 
ing, produc i ng mo r e early-type galaxie s as time passes 
Park et al.l mM IPark fc ChoiT l2(M iHwang fc Parkl 
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20081 ) . The speed of this process is an increasing function 
of the large-scale density. Namely, even though the di- 
rect physical processes affecting galaxy morphology are 
the gravitational and hydrodynamic interactions between 
neighboring galaxies, the large-scale background density 
appears to control galaxy morphology through its sta- 
tistical correlation with the frequency and strength of 
galaxy-galaxy interactions. Given the knowledge that 
the MDR in most region of the universe (note that Park 
et al. did not resolve the cluster regions) is the result 
of galaxy-galaxy interactions, one can naturally suspect 
the MRR and MDR in clusters are also due to the inter- 
actions between individual galaxies. It is the purpose of 
this paper to explore this possibility. 

2. OBSERVATIONAL DATA SET 

2.1. Sloan Digital Sky Survey Sample 

We use a spectroscopic sample of galaxies in the 
Sloan Digital Sky Survey (SPSS) D ata Release 6 
(DR6; lAdelman-McCarthv et all l2008fh The survey 
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Fig. 1. — Completeness of the spectroscopic sample of cluster 
galaxies as a function of r-band magnitude (upper panel) and clus- 
tercentric radius (lower panel). Filled circles are the spectroscopic 
completeness of the data complemented by NED, while open circles 
show the completeness of the original SDSS spectroscopic sample. 

produced five-band (ugriz) photometric data for 230 
million objects over 8,400 deg 2 , and optical spectro- 
scopic data more than one million objec ts of galax- 
ies, quasars, and stars over 6860 deg 2 dGunn et al.l 



les, quasars, ana stars o ver b »bu cleg tlljunn et al 
19981 120061 lUomoto et all 119991: ICastander et all 12001 



Blanton et alJ 120031: iFukugita et al.l 119961: iLupton et al 



2002t iHogg et all 1 20011: [Smith et all 120021 : llvezic et al 
12004 iTucker et al.ll2006l : (Pier et al.ll2 003f). E xtensive de- 
script ion o f SDSS data products i s given bv lYork et al.l 
(2000J) and lStoughton et al.l (|2002f ). 

The data is supplemented by several value-added 
galaxy catalogs (VAGCs) drawn from SDSS data. Photo- 
metric and structure para meters of galaxies are o btained 
from the SDSS pipeline (jStoughton et alj|2002j) . Com- 
plementary photometric parameters such as color gradi- 
ent, concentration index, and Petrosian ra d ius ar e taken 
from the DR4plus sample of IChoi et all (|2007| ). The 
spectroscopic para meters are obtained from MPA/JHU 
and N YU VAGCs (|Tremonti et all 12004 iBlanton et all 
l2005h . 

Completeness of the spectroscopic data in SDSS is poor 
for bright galaxies with m r < 14.5 because of the prob- 
lems of saturation and cross-talk in the spectrograph, 
and for the galaxies located in high density regions such 
as galaxy clusters due to the fiber collision (two spec- 
troscopic fibers cannot be placed closer than 55" on a 
given plate). Therefore, it is necessary to supplement 
the galaxy data to reduce the possible effects of the in- 
completeness problem. We search for the galaxies within 
ten times the virial radius of each galaxy cluster in the 
photometric catalog of the SDSS galaxies, and find their 
rcdshifts from the NASA Extragalactic Database (NED) 
to supplement our spectroscopic sample. Figure Q] shows 
the spectroscopic completeness of our galaxy sample as 
a function of apparent magnitude and of clustercentric 
distance. The open circles are the completeness of the 
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original SDSS sample, and the filled circles are that of 
our sample with additional redshifts. It shows that the 
spectroscopic completeness of our sample is higher than 
85% at all magnitudes and clustercentric radii. 

2.2. Cluster Sample and Galaxy Membership in 
Clusters 

We used the Abell catalog of galaxy clusters 
(|Abell et al.1 11989ft to identify cluster galaxies in our 
galaxy sample. Among the Abell clusters, we selected 
those that have known spectroscopic redshifts in the 
NED. We found 730 clusters located within the SDSS 
survey region. We adopted the position of cluster center 
in the NED, but replaced it with the X-ray determined 
position if it is available in the literature. 

In order to determine the membership of galaxies 
in a cluster, we used the "shifting gapper" method of 
iFadda et all (|1996t ) that was used for the study of k ine- 
matics of galaxy clusters (jHwang fc Ledl2007L I2008T ). In 
the radial velocity versus clustercentric distance space, 
the cluster member galaxies were selected by grouping 
galaxies with connection lengths of 950 km s _1 in the 
direction of the radial velocity and of 0.1 /i _1 Mpc in the 
direction of the clustercentric radius R. If the boundary 
was not reached out to R = 3.5/i~ 1 Mpc, we stopped the 
grouping at R = 3.5/i _1 Mpc. We iterated the procedure 
until the number of cluster members converges. From 
this procedure we obtained 200 Abell clusters that have 
more than or equal to 10 member galaxies. 

We computed a radius of r2oo,ci (usually called the 
virial radius) for each cluster where the mean overdensity 
drops to 200 times the criti cal density of the univ erse p c , 
using the formula given bv lCarlberg eTafl (fj~997h : 
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where a c \ is a velocity dispersion of a cluster and the 
Hubble parameter at z is H 2 (z ) = Hl[Q m (l + z) 3 + 
ti k (l + zf + Q A ] (|Peebleslll99"3]) . Q m , Q k , and Q A are 
the dimensionless density parameters. The velocity dis- 
persion was computed for each cluster from the redshift 
distribution of the cluster member galaxies as described 
in Appendix. 

In addition to the sample of cluster member galaxies 
obtained by the "shifting gapper" method above, we in- 
cluded the galaxies located at projected separations of 
Rmax < R < 10r2oo,ci to investigate the variation of 
galaxy properties over a wide range of clustercentric ra- 
dius. Rmax is the largest clustercentric distance of the 
cluster member galaxies determined above. These ad- 
ditional galaxies were constrained to have velocity dif- 
ference relative to the cluster's systematic velocity less 
than Aw = |w ga i — v sys \ — 1000 km s -1 . The final sample 



u sys\ 



consists of galaxies smoothly distributed from the clus- 
ter center to R = 10r2oo,ci for each cluster. Figure [5] 
shows the radial velocities of galaxies around eight clus- 
ters in our sample as a function of clustercentric distance 
of galaxies. 

We rejected the clusters that appeared to be inter- 
acting or merging, which was decided in the galaxy 
velocity versus clustercentric distance space. Dynami- 
cally young clusters having the brightest cluster galaxy 
(BCG) at large clustercentric distance (i?BCG > 0.18 hr 1 
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Fig. 2. — Radial velocity vs. clustercentric distance of galaxies. 
Filled circles and crosses indicate the early and late types, respec- 
tively, selected as galaxies associated with clusters, while open cir- 
cles indicate the galaxies not selected as associated members. The 
horizontal dashed lines indicate the systemic velocity of the clus- 
ters determined in this study. The vertical dotted lines indicate 
the radius r2oo,cl computed in this study. 

Mpc) were rejected too. We also eliminated the clus- 
ters for which survey coverages were not complete out 
to 10r2oo,ci- We finally obtained a sample of 93 relaxed 
Abell clusters and 34,420 associated galaxies for our anal- 
ysis. 

2.3. Physical Parameters of Galaxies 

The physical parameters of galaxies that we consider 
in this study are r-band absolute Petrosian magnitude 
(M r ), morphology, axis ratio, u — r color, equivalent 
width of Ha emission line, g — i color gradient, concen- 
tration index (c- m ), internal velocity dispersion (cr), and 
Petrosian radius in i-band. Here we give a brief descrip- 
tion of these parameters. 

The r-band absolute magnitude M r was computed us- 
ing the formula, 

M r = to,. - 51og[r(z)(l + z)] - 25 - K(z) - E(z), (2) 

where r{z) is the comoving distance at redshift z in unit 
of ft, Mpc, and the corresponding 51og/i term in M r will 
be omitted in this paper. K(z) is the if-correction, 
and E(z) is the luminosity evolution correction. Wc 
adopt a flat ACDM cosmology with density parameters 
Qa = 0.73 and fi m = 0.27. The rest-frame absolute 
magnitudes of individual galaxies are computed in fixed 
bandpasses , shifted to z = 0.1, using Galactic reddening 
correction ( Schlcgcl et al.)fl99 8) and if -corrections as de- 
scribed by iBlanton et al.l (120031). The evolution correc- 
tion given bv lTegmark et all (f200l . E(z) = 1.6(z-0.1), 
is also applied. 

Figure [3] shows the r-band absolute magnitudes of the 
cluster galaxies against their redshifts. We define the 
volume-limited samples of galaxies using the redshift and 
absolute magnitude conditions as follows: CI (—17.0 > 
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Fig. 3. — Sample definitions of our three volume-limited samples 
in the absolute magnitude vs. rcdshift space. The bottom curve 
corresponds to the apparent magnitude limit of m r = 17.77. 



M r > -19.0 and z < 0.0331), C2 (-19.0 > M r > -20.5 
and z < 0.0767), and C3 (-20.5 > M r > -22.5 and 
z < 0.1365). The apparent magnitude limit line (m r 
1 7.771 shown in Figure [3] is obtained using th e mean K- 
correction relation given by equation (2) of IChoi et al.l 
(|200l . 

Accurate morphology classification is critical in this 
work since the effects of interaction depend strongly on 
morphology of the target and neighbor galaxies. We 
first classify morphologi cal types of galax ies included in 
the DR4plus sample of IChoi et al.l (|2007f ) adopting the 
automated classification method given bv lPark fc Choi 
(2005). Galaxies are divided into early (elliptical and 
lenticular) and late (spiral and irregular) morphological 
types based on their locations in the u — r color ver- 
sus g — i color gradient space and also in the z-band 
concentration index space. The resulting morphologi- 
cal classification has completeness and reliability reach- 
ing 90%. The automatic classification scheme does not 
perform well when an early-type galaxy starts to overlap 
with another galaxy. This is because the scheme excludes 
galaxies with very low concentration from the early-type 
class and blended images often erroneously give low con- 
centration. Since we are investigating the effects of close 
galaxy-galaxy and galaxy-cluster interactions on galaxy 
properties, this problem in the automatic classification 
has to be remedied. We perform an additional visual 
check of the color images of galaxies to correct misclas- 
sifications by the automated scheme. In this procedure 
we changed the types of the blended or merging galax- 
ies, blue but elliptical-shaped galaxies, and dusty edge- 
on spirals. In addition, for the galaxies in DR6 that are 
not in the DR4plus sample, we visually classified their 
morphological types using the color images. 

The (it — r) color was computed using the extinc- 
tion and X-corrected model magnitude. The superscript 
0.1 means the rest-frame magnitude i^-corrected to the 
redshift of 0.1, and will subsequently be dropped. 

We adopt the values of (g — i) color, concentration 
index (cj n ), and Petrosian radi us Rp e t comput ed for 
the DR4plus sample of galaxies ()Choi et al. 2007). The 



(g — i) color gradient was defined by the color difference 
between the region with R < 0.5i?p c t and the annulus 
with 0.5i?p ct < R < Rpet, where i?p ct is the Petrosian 
radius estimated in z-band image. To account for the 
effect of flattening or inclination of galaxies, elliptical 
annuli were used to calculate the parameters. The (in- 
verse) concentration index is defined by i?5o / Rgo , where 
i?50 and Rqq are semimajor axis lengths of ellipses con- 
taining 50% and 90% of the Petrosian flux in the z-band 
image, respectively. 

The velocity disp ersion value of the g alaxy is adopted 
from NYU-VAGC (jBlanton et al.l 12003 ). The value of 
Ha equivalent width is taken from MPA/JHU-VAGC 
(|Tremonti et alj I2004D . which is computed using the 
straight integration over the fixed bandpass from the 
continuum-subtracted em ission line with the model of 
iBruzual & Charlotl (|2003h . 

In our analysis we often limit the late-type galaxy 
sample to galaxies with z-band isophotal axis ratio b/a 
greater than 0.6. This is to reduce the effects of internal 
extinction on our results. The absolute magnitude and 
color of late-type galaxies with b/a < 0.6 are very inaccu- 
rate (see Figs. 5 and 12 of lChoi et ai1l2007f K Therefore, 
including them in the analysis may introduce a large dis- 
persion in SF indicators such as luminosity, color, Ha 
equivalent width, and color gradient. 

2.4. Nearest Neighbor Galaxy in Clusters 

To account for the effects of the nearest neighbor 
galaxy in cluster environment, we determine the distance 
and the morphology of the nearest neighbor galaxy. 

We define the nearest neighbor galaxy of a target 
galaxy with absolute magnitude M r as the one which 
is located closest to the galaxy on the sky and is brighter 
than M r + AM r among those in our cluster galaxy sam- 
ple. We adopt AM r = 0.5. When we adopt galaxies 
fainter than the target galaxy by more than 0.5 mag as 
neighbors, our conclusions do not change but our statis- 
tics are worse since the number of target galaxies be- 
comes smaller (see ^3. II for more discussion). 

We do not use the velocity condition to determine the 
nearest neighbor galaxy because it is selected from the 
cluster galaxy sample to which the velocity condition is 
already applied. We obtain the nearest neighbor distance 
normalized by the virial radius of the nearest neighbor 
as follows. We first compute the small-scale density ex- 
perienced by a target galaxy attributed to its neighbor, 



Pn/p = 7„i„/(47rr^p/3), 



(3) 



where j n is the mass-to-light ratio of the neighbor galaxy, 
L n is the r-band luminosity of the neighbor, r p is the pro- 
jected separation of the neighbor from the target galaxy, 
and p is the mean density of the universe. We assume 
that 7(early)=27(late) at the same r-band luminosity, 
and that 7 is constant with galaxy luminosity for a given 
morphological type. The value of mean density of the 
universe, p = (0.0223 ±0.0005)(7L)_ 2 o(^" 1 Mpc)- 3 , was 
adopted, where ( 7^)-2o is the ma ss of a late-type galaxy 
with M r = -20 (jPark et al.ll2008ft . 

Then, we define the virial radius of a neighbor galaxy 
as the projected radius where the mean mass density p n 
within the sphere with radius of r p is 200 times the criti- 
cal density or 740 times the mean density of the universe, 
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Fig. 4. — The surface number density of early-type (upper 
panel) and late-type (middle panel) galaxies in the Abell cluster 
regions, and the corresponding early-type fraction (bottom panel) 
as a function of the clustercentric radius normalized to the cluster 
virial radius. The lines in the bottom panel are the best-fitting 
functions. 

namely, 

r vir = (3 7 L/4 7 r/200 jOc ) 1 / 3 . (4) 

Since we adopt O m = 0.27, 200p c = 200 / o/fi m = 740p. 
This is almost equal to the virialized density p v irial = 
l8ir 2 /n m (H Q t ) 2 p = 7665 in the case of our ACDM 
univer se (jGott fe ReesHl975f) . This is what iPark et all 
(2008) used to define the virial radius. According to 
our formula the virial radii of galaxies with M r = 
-19.5,-20.0, and -20.5 are 260, 300, and 350 hr 1 kpc 
for early types, and 210, 240, and 280 hr 1 kpc for late 
types, respectively. 

3. RESULTS 
3.1. Morphology-Environment Relation 

Figure 2] shows the surface number density of galax- 
ies and the fraction of early-type galaxies as a function 
of projected clustercentric radius R normalized by the 
cluster virial radius 7"200,ci- They are shown for three 
luminosity ranges. The uncertainties of the fraction rep- 
resent 68% (la) confidence intervals that are determined 
from the numerical bootstrap procedure. 

It shows the well-known MRR or MDR of galaxies 
in clusters. In all three luminosity ranges the early- 
type fraction Je is almost constant at large distances, 
and starts to increase inwards at the critical region with 
R rj 1 ~ 3r2oo,clj which corresponds to the cluster in- 
fall region. The transition radius seems depending on 
galaxy luminosity. For relatively brighter galaxies with 
—20.5 > M r > —22.5 (the subsample C3) the transi- 
tion occurs at ~ l-5r2oo,ci, but for fainter galaxies with 



4.0 L""l I I I 




1.0 r .... I I I I " 

-21 -20 -19 -18 -17 
M r (mag) 

Fig. 5. — The characteristic clustercentric radius R CI as a func- 
tion of r-band absolute magnitude. R cr is the scale where the early- 
type fraction Je starts to rise significantly above the field value. 
Five luminosity subsamples are used: —17 > M r > —18, —18 > 
M r > -19,-19 > M T > -19.75,-19.75 > Mr > -20.5,-20.5 > 
M r > -22.5. The best-fit linear line is fl C r/»'200,cl = (0.62 ± 
0.10)M r + (14.32 ± 1.95). 

— 17.0 > M r > —19.0 (the subsample CI) it occurs at 
~ 3r2oo,ci- The slope of increase is larger for less lumi- 
nous galaxies at small R. Since we are not distinguishing 
between ellipticals and lcnticulars, we can only see their 
sum is rising toward the cluster center. Fainter galaxies 
seem more vulnerable to the cluster influence. This can 
be considered as evidence for the direct interaction be- 
tween clusters and infalling galaxies, thinking that less 
massive galaxies change their morphology at farther dis- 
tances from clusters relative to massive ones. However, it 
can be also due to the mass segregation within clusters 
that massive cluster member galaxies hardly overshoot 
beyond the cluster virial radius while the less massive 
ones overshoot out to a few times the virial radius. Mas- 
sive galaxies outside the virial radius are likely to be the 
infalling ones rather than bound cluster members. 

Figure 3] indicates that the early- type fraction is nearly 
the same near the cluster center regardless of galaxy lu- 
minosity. The central value of Je = 0.8 ~ 0.85 actually 
means that the fraction is close to 1.0 at the cl uster center 
if the projection effect is taken into account (|Ann et al.l 
2008). At radii larger than about 4r2oo.ci 7 /b approaches 
the field fraction that depends on the luminosity range 
and becomes inse nsitive to the large-scale environment 
(|Park et al.ll200l . 

To inspect this characteristic scale in more details we 
fit Je as a function of i?/r2oo,ci using the function 

Ie = (0.9 - J Q )cM-R/Ro) + fo, (5) 

where Rq is a free parameter and /o is an average value 
of Je at -R/ r 200.ci > 4.0. We define the characteristic 
radius R cr as the scale where Je becomes 10% larger 
than its field value at fixed luminosity. Figure [5] shows 
an almost linear relation between the characteristic ra- 
dius and the absolute magnitude obtained from five lu- 
minosity subsets. The best-fit linear line is -Rcr/^oo ci = 
(0.62 ± 0.10)M r + (14.32 ± 1.95). 

We now investigate the dependence of galaxy morphol- 
ogy on both clustercentric radius and nearest neighbor 
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Fig. 6. — Morphology-environment relation when the nearest galaxy is (left) an early-type and (right) a late-type. Absolute magnitude 
of galaxies is fixed to a narrow range of —19.5 > M r > —20.0. Contours show constant early-type galaxy fraction fg. The points with 
error bars above the x-axes denote the average virial radius of the BCGs. 



distance. We select a volume-limited sample of target 
galaxies associated with clusters whose absolute magni- 
tudes are in a narrow range of —19.5 > M r > —20.0 and 
rcdshifts are less than 0.0767. Their nearest neighbors 
are found among all galaxies with M r < —19.0. Then 
r v j r of the nearest neighbors are calculated using their 
luminosity and morphology. 

Dots in Figure [5] show the distribution of early-type 
(red points) and late- type (blue points) galaxies in the 
projected clustercentric radius R and projected nearest 
neighbor distance R n space. A spline kernel is used to 
obtain the smooth distributions of the median Je in each 
location of the two panels. Contours with different col- 
ors mark constant early-type fractions. The left panel of 
Figure [5] contains the target galaxies having early-type 
neighbors. It shows that at R > 2r2oo,ci ah contours are 
nearly horizontal. This means that outside the cluster 
virial radius galaxy morphology is determined solely by 
the nearest neighbor distance and morphology. When the 
nearest neighbor is an early type, fs monotonically in- 
creases as R„, decreases. But if the neighbor is a late type, 
it first increases, reaches a maximum at R n ~ r v i r , n /3, 
and then decreases as R n decreases. The bifurcation oc- 
curs at Rn ~ fVi^n. This is the morpholo gy-neighbor 
enviro nment relation that was discovered by 1 Park et alJ 
(2008) in the general large-scale background density en- 
vironment. It is interesting to see that the effects of the 
nearest neighbor's distance and morphology are the dom- 
inant factors of galaxy morphology transformation right 
down to the cluster infall regions. 

The situation abruptly changes at the critical clus- 
tercentric radii of i? cr = 1 ~ 3r2oo,ci- Within this ra- 
dius contours arc nearly vertical when R n > O.lrvi^n, 
but are nearly horizontal when R n < 0.05r v j r!n . When 
Rn ^ 0.1rvir,n> Ie monotonically increase as R decreases 
and depends almost entirely on R regardless of the mor- 
phological type of the nearest neighbor. Existence of this 
sudden transition near the cluster virial radius suggests 
that the morphology transformation of cluster galaxies 
is due to the interactions between galaxies and clusters. 



Both gravitational and hydrodynamic processes can pro- 
duce the discontinuity in Je- It may be the hot cluster 
gas confined within ~ »"200,ci or gravitational tidal force 
of cluster acting on galaxies trapped within the virial 
radius that causes the transformation. The monotonic 
increase of Je at smaller R can arise by the increase 
of the effects of hot cluster gas whose pressure mono- 
tonically increases toward the cluster center and/or by 
interactions with the cluster potential for galaxies with 
smaller orbital radii that have the shorter crossing time. 

On the other hand, Figure [6] clearly indicates that 
the clustercentric radius is not the only environmental 
parameter determining the galaxy morphology but the 
nearest neighbor does a critical role when the neigh- 
bor separation is less than about 0.05r v ; rin . The local 
galaxy number density, to which the MDR is often at- 
tributed, cannot be responsible for the increase of /e 
inwards cluster because the density rises smoothly as R 
decreases and does not have a characteristic break at R CI 
(see Fig. [4j. After all, it is the morphology-clustercentric 
radius-neighbor environment relation instead of the sim- 
ple MDR or MRR that determines the morphology of 
cluster galaxies. Here neighbor environment includes 
both neighbor distance and neighbor morphology. 

We emphasize that we are not using the nearest neigh- 
bor distance as a measure of local galaxy number density. 
It is a measure of the influence of the nearest neighbor 
galaxy itself irrespectively of other galaxies. However, 
there exists statistical correlation between R n and the 
local galaxy number density. If one uses a measure of 
local galaxy number density, one will still find some cor- 
relation of the measure with galaxy morphology within 
the cluster virial radius. But it is actually the galaxy- 
galaxy interaction that causes the correlation. This is 
because fixing R essentially fixes the local galaxy num- 
ber density and in Figure [5] Je still changes at a fixed 
R. 

In the left panel of Figure [6] galaxies have early-type 
neighbors. It is expected that both cluster and neighbor 
galaxy impose morphology transformation toward early 
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type on galaxies through hydrodynamic or gravitational 
effects. Dominance of the role between the two depends 
on R and R n . At Rta r 20 oxi the neighbor galaxy starts to 
give a major impact on galaxy morphology at the separa- 
tions R n < 0.5r V i r n = 100 ~ 150/i _1 kpc (note the differ- 
ence of Je between the two panels). But at R « 0.2r 2 oo, c i 
it happens at R n < 0.05r virin « 10 ~ 15/i _1 kpc. This 
seems reasonable since near the cluster center the tidal 
effects by the cluster potential are stronger, the cluster 
gas has higher pressure, and thus the cluster gives more 
direct impact on galaxy properties. 

One can note from a comparison between two panels 
of Figure [6] that directions and levels of contours change 
completely depending on the morphology of the nearest 
neighbor galaxy at any clustercentric radius when the 
pair separation is small (i.e. R n < 0.1r v ; r n ). The right 
panel shows that a cluster galaxy tends to become a late 
type when it has a late-type neighbor within the sepa- 
ration R n w 0.1r v i rjn even when the pair is well within 
the cluster virial radius. For pairs with R n « 0.05r v ; r n 
Je ~ 0.8 at the clustercentric radius of R = 0.5r 20u , c i 
when the neighbor is an early type. But Je is as low as 
0.35 for pairs with the same R n located at the same R 
when the neighbor is a late type. It is a clear demonstra- 
tion that the nearest neighbor has a dominant control 
over galaxy morphology transformation in clusters when 
the distance to the neighbor is less than about one tenth 
of the virial radius of the neighbor galaxy. 

Figure [6] shows that there are not many such short 
separation pairs at a particular moment within the clus- 
ter virial radius and one might think that the majority 
of cluster galaxies are not affected by neighbors. How- 
ever, it can be seen in Figure [5] that at R < 0.5r 2 oo, c i 
the mean separation between galaxies becomes so small 
that the virialized regions associated with galaxies are 
all overlapping with one another (not that the upper 
edge of the galaxy distribution becomes smaller than 
Rn/r V ir,n = 1 at i?/r 2 oo,ci < 0.5 in Fig. [5]). Suppose the 
RMS orbital velocity of galaxies within r2oo,ci is {v 2 ) 1 / 2 . 
If r 20 o,ci = 1.5 h- 1 Mpc and (v 2 ) 1 / 2 = 700 km s -1 , 
the crossing time across r 2 oo, c i is about 3 x 10 9 yrs. On 
the other hand, the RMS relative velocity between two 
random cluster galaxies is ^{v 2 ) 1 / 2 . If the virial ra- 
dius of the galaxies is r v i r = 300 hr 1 kpc, their cross- 
ing time is 4 x 10 8 yrs, an order of magnitude shorter 
than the crossing time across the cluster. This means in 
this figure that the cluster galaxies at R < r 2 oo,ci have 
migrated vertically (interaction with neighbor galaxies) 
many times during the time they make one travel hor- 
izontally (interaction with the cluster). In particular, 
when R < 0.5r 2 oo, c i, a galaxy starts to interact with an- 
other neighbor galaxy as it passes one neighbor galaxy 
because the virial radii of galaxies all overlap with one 
another there. Then the nearest neighbors may have 
left significant cumulative effects on morphology of clus- 
ter galaxies. The amount of the cumulative effects in- 
creases monotonically as the clustercentric radius de- 
creases, which can produce the trend seen in Figure [5] 

It is expected that the late-type galaxy sample has a 
higher fraction of interlopers in the cluster environment. 
The pairs involving late- type target or late-type neighbor 
are more likely to be false pairs, and more widely sepa- 
rated in real space. If this effect is corrected in Figure 
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Fig. 7. — Absolute magnitude of galaxies brighter than M r = 
— 19 in the Abell cluster regions as a function of clustercentric 
radius. The upper panel shows the early types, and the lower panel 
shows the late types. Lines are the median magnitudes. Crosses 
are the BCGs, and they are not used in calculating the median 
values. Late types with axis ratio of b/a < 0.6 were eliminated. 
The vertical dashed lines indicate the cluster virial radius T2oo cl • 



O the contours will be shifted upward. Furthermore, if 
the fraction of interlopers varies as a function of neigh- 
bor separation or clustercentric radius, the contours will 
be both shifted and stretched. However, the fact that 
there are discontinuities in ]e shown in Figure [6] (near 
the cluster virial radius R rs t" 2 oo jC i and the galaxy merger 
scale R n rts 0.05r v i rjn ) and also the fact that these tran- 
sitions occur at physically meaningful scales, cannot be 
explained if the fraction of interlopers is seriously high. 
This suggests that the morphology-environment relation 
shown by the contours in Figure[6]is not seriously affected 
by the interlopers. 

Though we rejected the clusters that appear to be in- 
teracting or merging, there are some clusters (~15%) 
whose projected distance to the nearest cluster is less 
than 10r 2 oo.ci and relative velocity is less than 1000 km 
s _1 . Therefore, the galaxies associated with these clus- 
ters have been counted more than once in our analysis. 
For example, among 3874 galaxies used in Figure [6j 874 
galaxies are those counted more than once. When we re- 
ject these galaxies, our results remain almost the same. 
In addition, our results do not depend sensitively on our 
definition of the nearest neighbor. For example, when we 
varied AM r from 0.5 to .75, our results remain qualita- 
tively the same (see also iPark et alj[2008t iPark fc Choi 
120091) . 

The mean virial radius of the BCGs, which is typically 
about 0.5r 2 oo,ci, and its la range are shown above the 
x-axes of Figure [H We do not find any characteristic fea- 
ture across the virial radius of the BCGs. This suggests 
that the BCGs are not directly participating in trans- 
forming the morphology of the cluster galaxies. When 
the positions of BCGs are used as the centers of clus- 
ters, the contours in Figure [5] are basically unchanged at 
R > 0.2r 20 o,ci- But the slight drop of Je at R < 0.1r 2 oo, c i 
seen for the E/S0 neighbor case, now disappears. 
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Fig. 8. — Median absolute magnitude contours in the projected pair separation R n /i" v i r>n vs. the clustercentric distance R/f"2oo,cl f° r 
the galaxies with z < 0.0767 and M r < —19.5. Late types with axis ratio of b/a< 0.6 were eliminated. Four cases are given; the early-type 
target galaxies having an early- type neighbor (E-e), the early- type target galaxies having a late- type neighbor (E-l), the late- type target 
galaxies having an early-type neighbor (L-e), and the late- type target galaxies having a late- type neighbor (L-l). 



3.2. Luminosity- Environment Relation 

Figure [7] shows the r-band absolute magnitudes of the 
galaxies in a volume-limited sample with z < 0.0767 and 
M r < —19.0. The lines with error bars are the median 
values as a function of R. The BCGs, marked as crosses, 
are not used in calculating the median values. Early-type 
galaxies become fainter as much as ~ 0.2 magnitude as R 
decreases from 10r2oo, c i to 0.2r2oo,cl- The luminosity of 
early types rises again as R becomes smaller than about 
0.2r2oo.ci- But the increase of luminosity near the cluster 
center seems to depend on how much the BCGs are off- 
centered. Such a variation of the median luminosity is 
not detected for late types. 

The dependence of luminosity on the environment can 
be better understood in the two dimensional study. Con- 
tours in Figure [5] delineate the constant median M r lo- 
cations in the two-dimensional environmental parameter 
space of R and R n . Four panels distinguish among four 
different combinations of target and neighbor morphol- 
ogy. At each location of the R-R n space the median value 
of M r is found using the galaxies included within a spline 
kernel with a fixed size. The BCGs are not included in 
the smoothing. 

It can be clearly seen that galaxies with larger R n are 
brighter when R is greater than about 0.2r2oo,ci- Namely, 
the bright galaxies are more isolated from their influen- 
tial neighbors than the relatively faint galaxies are. In 
the region far outside the cluster virial radius the abso- 
lute magnitude of galaxies is nearly constant in R (i.e. 
the contours are nearly horizontal). Luminosity of galax- 
ies with R n > r v i rjn depends very sensitively on R n . But 
galaxies with R n < r v ; rin show little dependence of M r 



on R or Rn. The same phenomenon has been found by 
iPark et all (|2008f) for galaxy pairs in the general back- 
ground density regions. 

The reason that galaxies are brighter when they are 
located outside the nearest neighbor's virial radius was 
interpreted by Park et al. as due to luminosity trans- 
formation through mergers. After a galaxy merges with 
its closest neighbor, the second nearest neighbor will be- 
come the new nearest neighbor of the merger product. 
Since the merger product is typically located at a sepa- 
ration from the originally second nearest neighbor larger 
than the virial radius of the neighbor, a recently merged 
galaxy will appear in general at the upper edge of the 
galaxy distribution in Figure [8j As one moves toward 
the cluster center, the average galaxy number density 
increases and the mean distance between galaxies de- 
creases. Correspondingly, the distance to the new near- 
est neighbor of a recently merged galaxy also decreases 
statistically for decreasing clustercentric radius, and the 
location of the recent merger products follows the upper 
edge of the galaxy distribution in Figure [8J 

The upper left panel of Figure [5J showing the early- 
type target galaxies having early-type neighbors (the E- 
e case), shows that the maximum luminosity a galaxy 
can attain slowly decreases as R decreases until R be- 
comes about 0.2r2oo,ci- The median M r of the most sep- 
arated galaxies is about —21.3 at R w 10r2oo,ci, but de- 
creases to about —20.0 at R w 0.2r2oo,ci, a more than 
one magnitude drop. This can be explained if the lumi- 
nosity transformation process by galaxy-galaxy mergers 
becomes less efficient as R decreases. This is reasonable 
because the relative velocity between neighboring galax- 
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Fig. 9. — Dependence of star formation activity parameters of our target galaxies in the samples of CI, C2, and C3 on the clustercentric 
distance: (left) u — r, (middle) W(Ha)+l(A), and (right) A(g — i). Median curves are drawn for the cases of early types (solid line) and 
of late types (dotted line). Late types with axis ratio of b/a < 0.6 were eliminated. 



ies will be higher and the chance for a pair to merge will 
decrease toward the cluster center. 

A similar independent trend in luminosity is also seen 
for late-type galaxies as shown in the right panels of Fig- 
ure [8] even though the trend is much weaker. Luminos- 
ity of late-type galaxies is on average lower than that 
of early types. Interestingly, the late-type galaxies hav- 
ing an early-type neighbor (the L-e galaxies) are slightly 
brighter than the early-type galaxies at R < 0.1r2oo,ci- It 
could be that, as R decreases, faint late types are trans- 
formed to early types and only relatively bright late types 
survive near the cluster center. 

Figure [5] also shows that galaxies located near R = 
T200, ci and having a close neighbor are on average the 
faintest. Namely, even though the most isolated galax- 
ies become faintest at R w 0.2r2oo,ci, those with a 
close neighbor (R n w 0.1r V i r , n ) become faintest at R « 
l-Or 2 oo,ci- 

3.3. Star Formation Activity Parameters 

The left column of Figure [9] shows the u — r color of 
galaxies divided into two morphology and three magni- 
tude bins as a function of R. It can be seen that the 
galaxy color becomes redder as R decreases for both early 
and late types. This clustercentric radius dependence of 
color is stronger for relatively fainter galaxies (i.e. the CI 
sample). In the case of the galaxies in CI the dependence 
seems to occur at ~ 3r2oo,ci- 

Dots in Figure [10] are the intermediate luminosity 
galaxies with —19.5 > M r > —20.5 in the two- 
dimensional environmental parameter space. Color of 
early-type galaxies is nearly constant everywhere. There 
is a slight tendency that early-type galaxies become red- 
der as they approach an early-type neighbor. 



The constant color contours of late-type galaxies are 
nearly horizontal at R > R cl - = 1 ~ 3r2oo,ci, which means 
galaxy color is independent of R and depends only on the 
neighbor separation R n . Inside the critical radius i? cr the 
color of late-type galaxies still depends mainly on neigh- 
bor separation and morphology. This is evidence that hy- 
drodynamic interactions between individual galaxies are 
effective even within the virialized region of clusters. In 
the case of the L-l galaxies there is clustercentric radius 
dependence of color when R n > r v i rjn /3. It is likely that 
this is because the neighbors other than the nearest one 
tend to be early types as the system moves toward the 
cluster center. If the cluster hot gas had a direct impact 
on late-type galaxies' color, both L-e and L-l galaxies 
would show a similar color dependence on R. 

The equivalent width of the hydrogen Ha line is a mea- 
sure of SFA. Figure [9] shows the median W(Ha) as a 
function of R for six different cases distinguished by mor- 
phology and luminosity of target galaxies. The SFA of 
late-type galaxies decreases as R decreases. The depen- 
dence again seems to start at ~ 3r 2 oo,ci for the galaxies 
in the CI sample. It is important to note that both 
u — r color and W(Ha) show only mild dependence on 
the clustercentric radius i?, and it is actually the neigh- 
boring galaxies that mainly determine the color and SFA 
of galaxies within clusters as we can see in Figure [TO] 

The second column of Figure [TO] shows W(Ha) in two- 
dimensional environmental parameter space. One can 
notice that the constant W(Ha) contours look very sim- 
ilar to those of u — r color. Namely, the SFA of the L-e 
galaxies is mainly controlled by the neighbor galaxies. 
On the other hand, the SFA of L-l galaxies shows de- 
pendence on R when the nearest neighbor separation is 
not too small (R n > 0.3r V i rjn ). This is again probably be- 
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Fig. 10. — Dependence of u — r color, equivalent width of the Ha line, g — i color gradient of galaxies with —19.5 > Mr > —20.5 on 
the pair separation R n and the clustercentric distance R. In each column, galaxies are divided into four cases, the E-e, E-l, L-e, and L-l 
galaxies. Dots are galaxies belonging to each subset. At each location of the i?n/f v ir.n-^/ r 200,cl space the median value of the physical 
parameter is found from those of galaxies within a certain distance from the location. Curves are the constant-parameter contours. Late 
types with axis ratio of b/a < 0.6 were eliminated. 



cause the neighbors other than the nearest one tend to be 
early types as R decreases. The two cases demonstrate 
that the SFA of the late-type cluster galaxies changes 
differently depending on the morphology of the neighbor 
galaxy even well inside the cluster virial radius. If the 
hot cluster gas could directly quench the SFA of late- type 
galaxies, both L-l and L-e galaxies would show the same 
i?-dependence of SFA regardless of the morphology of the 
nearest neighbor. Since this is not the case, one can con- 
clude that it is after all the galaxy-galaxy hydrodynamic 
interaction that gives the biggest impacts on the color 
and SFA of cluster galaxies and that, contrary to intu- 
ition, the hot cluster gas is not the main actor quenching 
the star formation in late-type cluster galaxies at least 
down to R « 0.1r2oo,ci- We have repeated our analyses 
using the fainter galaxies with —17.5 > M r > —19.0, and 



arrived at the same conclusion. In the case of early-type 
galaxies the W(Ha) parameter is very small and hardly 
changes. 

Figure [5] shows that the color gradient of galaxies de- 
pends only weakly on R. In all cases the central region of 
galaxies becomes slightly bluer relative to the outskirts 
as R decreases. The dependence is stronger for fainter 
galaxies. 

In the right column of Figure flOl contours represent the 
distribution of the median A(g — i) at each location of 
the R-R n space. Unlike the u — r color and W(Ha), the 
color gradient A(g — i) of late- type galaxies depends on 
both R and R n inside r2oo,ci as can be seen from the slant 
contours. We interpret this phenomenon as a result of 
galaxy-galaxy interaction. When a late-type galaxy ap- 
proaches a neighbor closer than the virial radius of the 
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neighbor, its color gradient always increases (center be- 
comes relatively bluer) regardl e ss of n eighbor's morphol- 
ogy as shown bv lPark fe; Ch oi (2009, see Figs. 6 and 7). 
The color gradient of a late-type galaxy can increase as 
it moves toward the cluster center because it becomes 
more likely to be affected by neighbor galaxies. Since 
there is no neighbor morphology dependence of A(g — i) 
when a late type approaches a neighbor, the A(g — i) 
contours are quite similar for L-e and L-l galaxies. It is 
also possible that the center of late-type cluster galaxies 
becomes bluer at small er R because of the tidal effects 
of the cluster potential (Merritt 1984). 

3.4. Structure Parameters 

The left column of Figure [TTJ shows that concentra- 
tion of intermediate and low luminosity late-type galax- 
ies increases (ci„ decreases) as R decreases below R CI . 
But concentration of early types and high luminosity 
late types is nearly independent of R. As argued by 
iPark fc Choil d2009fl . this can be attributed to less com- 
pact internal structure of late types which are more vul- 
nerable to tidal effects than early types. 

The first column of Figure [T2] shows the median Cj„ 
contours for the intermediate luminosity galaxies with 
— 19.5 > M r > —20.5. Again ct n of early types hardly 
changes as R or R n vary. On the other hand, Cj„ of late 
types starts to decrease (galaxies become compacter) at 
R < R CI = 1 ~ 3r2oo,ci- It is interesting to see that the 
structural parameter Ci n shows an abrupt change at the 
characteristic scale. A similar phenomenon can be found 
also for the central velocity dispersion. 

The fact that the i?-dependence of Ci n exists only when 
R < i? cr , gives an important clue for understanding the 
structural evolution of cluster galaxies. The local galaxy 
number density cannot be the reason for the discontin- 



uous i?-dependence of Ci„ because the local density is a 
smooth function of R. The discontinuity can appear due 
to repeated gravitational interactions of cluster member 
galaxies with cluster potential or with other galaxies as 
they make trapped orbital motions within the cluster 
virial radius. It might also seem possible to explain the 
discontinuity by the hydrodynamic effects of hot cluster 
gas. For example, the cold gas in the outer part of late- 
type galaxies can evaporate or be stripped when they fall 
into the hot cluster gas clump or encounter the hot halo 
gas of their neighbor galaxies. Then their outer part will 
become redder and dimmer. This may cause the color 
u — r and color gradient A(g — i) increase and inverse 
concentration index c„ decreases as R or R n decreases. 
However, it is difficult to explain why the central stel- 
lar velocity dispersion of late types should also increase 
shown in the second column of Figure [T^] by quenching 
the SFA in the outer part. 

Figure [TTJ (panels d, e, and f ) shows that the central 
velocity dispersion of galaxies in general increases as R 
decreases except for low luminosity galaxies. The trend is 
clear for intermediate luminosity galaxies, in particular. 

The second column of Figure [T2] shows constant a con- 
tours of intermediate luminosity (—19.5 > M r > —20.5) 
galaxies divided into four different target-neighbor mor- 
phology cases. The central stellar velocity dispersion a 
shows a characteristic behavior that it depends only on 
Rn at R > i? cr but starts to depend only on R when 
R < R CI . Would the increase of a within the cluster 
virial radius be due to the tidal interaction of galaxies 
with other galaxies or with the cluster potential? One 
can find a clue from Figure [TTJ which shows that the 
i?-dependence of Ci n and a is largest for intermediate 
luminosity galaxies. This may be because the interme- 
diate luminosity galaxies have suffered from the interac- 
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Fig. 12. — Same as Fig. 1101 but for (left) Ci n , (middle) a (km s 1 ), and (right) Rp c t4 (h 1 kpc). Late-types with axis ratio of b/a < O.i 
were included. 



tion with cluster potential more strongly than brighter 
or fainter ones. Even though galaxy-galaxy encounters 
in general increase a of lat e-type galaxies (see L-e cases 
at R > r 20 o,ci and Fig. 8 of iPark fc C hoi 2009), it is not 
true at R < r 2 oo,ci as can be seen in Figure [TJ] The very 
bright massive galaxies at 0.1 < i?/ r 200,ci ^5 1-0 may be 
the galaxies falling into the cluster for the first time, and 
those which fell in long time ago are now sitting near 
the cluster center. In both cases the effects of the cluster 
potential will be relatively small. On the other hand, the 
small faint galaxies are likely to have large orbital radii, 
and can cross the cluster much fewer times than the mas- 
sive ones. This is supported by the fact that i? cr is larger 
for fainter galaxies (see Fig. [5]for example). Therefore, it 
is expected that the intermediate luminosity galaxies are 
most susceptible to change in structure and kinematics 
through interactions with the cluster potential. 

Our study of Ci n and a of cluster galaxies makes it 



clear that late-type galaxies become more centrally con- 
centrated and have higher central velocity dispersion as 
they approach the cluster center. Our results suggest 
that late-type galaxies become earlier in SFA through 
hydrodynamic interactions between galaxies and also in 
morphology and kinematics through tidal interactions as 
they approach the cluster center. Paucity of late-type 
galaxies near cluster center is a result of combined ef- 
fects of gravitational and hydrodynamic interaction with 
the cluster and the nearest neighbor galaxies. 

We use the Petrosian radius in the ?-band image as 
a measure of galaxy size. Figure [TT] (panels g, h, and 
i) shows that the intermediate luminosity galaxies show 
very slight decrease of i?p c t as R decrease. The third 
column of Figure [12] shows again that there is a weak 
clustercentric radius dependence of i?p e t at R < R cr for 
late-type galaxies. The decrease of i?p e t does not seem to 
be due to the blending of galaxy images in high density 



Galaxy Interactions in Clusters 



13 



regions like clusters because it suddenly occurs at the 
physically meaningful scale (i? C r)- The degree of blending 
of galaxy images should be a smooth function of R. 

The size of late-type galaxies does not change signif- 
icantly by interactions with other galaxies as shown in 
Figure 8 of iPark et all (|2008h . Therefore, the decrease of 
Rp e t of late types at R < R cl - shown in Figures [TT1 and [T2l 
should be attributed to the cluster. Late-type galaxies 
can appear smaller when their outer part becomes dim- 
mer as they orbit within a cluster and get quenched in 
star formation in disks. 

4. DISCUSSION 

Many mechanisms have been proposed to explain the 
MRR or MDR in clusters. They are divided into two 
categories : mechanisms re lying on gravitational an d hy- 
drodynamic processes (see iBoselli fc Gav azzi 2006| for a 
review). The mechanisms based on gravitational inter- 
actions include galaxy-galaxy tidal interaction, galaxy- 
cluster potential tidal interaction, and interaction with 
the general tidal background (harassment). Those based 
on hydrodynamic interactions comprise ram pressure 
stripping, viscous stripping, or thermal evaporation of 
cold gas of late-type galaxies by hot cluster gas, and re- 
moval of hot halo gas reservoir causing stopping of gas 
supply and quenching of star formation (strangulation). 
We will examine each of these mechanisms against our 
findings. 

1. Galaxy-galaxy tidal interaction 

Tidal interaction between galaxies can be efficient at 
removing or consuming the cold gas in galaxies, and tend 
to transform late t ypes i nto early types (ISpitzer fc Baade j 
1951 1 iRichstond Il976t iFarouki fc Shapiro! 119811: llckel 
1985 : [MerrittJll983h . However, this mechanism is usually 
excluded as a process responsible for the MRR or MDR 
because the orbital velocity of cluster galaxies are very 
high and the tidal energy deposit during the short en- 
coun ters is too small to significantly affect galaxy struc- 
ture (lMerrittlll984l : lBvrd fc Valtonenlll990f ). 

At R > R CT , the structural parameter of late- type 
galaxies depend only on the environment determined by 
the neighbor galaxies. But R < R cr they strongly de- 
pend on R, but only weakly or negligibly on R n (see Fig. 
T2"|) , This tells that the galaxy-galaxy tidal interactions 
do not instantly change the structure of cluster member 
galaxies much. The parameters Ci n and i?p e t do show a 
weak dependence on R n for L-e galaxies within R = R cr . 
But this may have also been resulted by hydrodynamic 
interactions with their early-type neighbor. Even though 
we exclude the significance of the tidal effects of galaxy- 
galaxy interactions on the internal structure and kine- 
matics of cluster galaxies, it should be emphasized that 
the hydrodynamic effects of galaxy-galaxy interactions 
can give significant impacts on their SFA. 

2. Galaxy-cluster potential interaction 

Individual galaxies can s uffer from the tidal force 
of the overall cluster mass dMerrittl [l98l iMilled 1198 
Bvrd fc ValtonCTl fl990l IMoss fe Whittle! l2000t IG^edl 
2003) . This mechanism is usually excluded in explaining 
the MRR or MDR because the time scale required for 
morphology transformation, being several cluster cross- 
ing time, is rather too long. The tidal force will be 
stronger toward the cluster center, and the strong de- 
pendence of Ci n and a on R within R = R CI can eas- 



ily explained by this mechanism. The internal structure 
of the intermediate mass galaxies, in particular, seems 
to have been significantly affected by this mechanism. 
But the tidal interaction with cluster potential may not 
be sufficient for morphology transformation of the very 
bright or faint cluster galaxies. 

3. Harassment 

Galaxies can be also perturbed by the tidal force back- 
ground from nu merous distant encounters wit h other 
cluster members (jMoore et al.lll996l I1998L |1999( ). Even 
though individual encounters are too short for the tidal 
force to activate structural changes, galaxy structure 
may significantly change through such numerous weak 
encounters. The sudden transition of the structure pa- 
rameters can be explained by repeated passage of mem- 
ber galaxies within the cluster. So both interaction with 
cluster potential and harassment may account for the R- 
dependence of morphology and structure parameters of 
intermediate mass galaxies, in particular, whose cluster 
crossing time is reasonably short but the path length is 
still long. 

4. Galaxy-galaxy hydrodynamic interaction 

The effects of hydrodynamic interaction between 
galaxies on cluster galaxy properties have been found 
important in this work. The galaxy-galaxy interaction 
turns out to be the major mechanism for quenching SFA 
of late-type galaxies within clusters. As can be seen 
in the scatter plot of galaxies in the R-R n plane, dis- 
tance to the neighbors monotonically decreases as R de- 
creases and galaxies become always located within neigh- 
bor galaxy's virial radius when R < i?200,ci/3. In this sit- 
uation galaxies are constantly undergoing hydrodynamic 
effects and mass exchange with neighbors. We argued 
above that the vertical contours for Ci n and a can be 
explained by the gravitational tidal force from cluster 
itself. However, the nearly horizontal contours for the 
u — r and W(Ha) parameters can be only explained by 
hydrodynamic effects of neighbor galaxies. It should be 
emphasized that the high-speed orbital motion of cluster 
galaxies weakens the importance of galaxy-galaxy tidal 
interactions, but can make the galaxy-galaxy hydrody- 
namic interaction very important. 

It was also found that bright galaxies are more isolated 
from influential neighbors than relatively fainter ones in 
clusters, which has been interpreted as due to merger- 
driven luminosity transformation. The fact that early- 
type galaxies show such trend more strongly is consis- 
tent with this interpretation. Relative paucity of bright 
galaxies near cluster center (except for the BCGs) seems 
to indicate that the frequency of galaxy merger decreases 
as the clustercentric radius decreases. 

5. Ram Pressure Stripping 

Many rich clusters like the Abell clusters we are an- 
alyzing, are holding hot X-ray emitting gas. A late- 
type galaxy falling into such a hot gas tank can strip 
off its cold gas an d terminate SFA (|Gunn fc Gott!H972l : 
IQuilis et alB oOO). The hot gas is trapped roughly within 
the virial radius of the cluster. Therefore, any mechanism 
relying on hydrodynamic interaction with the hot clus- 
ter gas will have a characteristic scale of onset of cluster 
influence near the cluster virial radius. 

A major problem with the ram pressure stripping and 
other hydrodynamic processes in accounting for the mor- 
phology segregation, is that these mechanisms can only 
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Fig. 13. — Same as Fig. [6] but for (a and b) X-ray bright clusters (Lx > 2 X 10 44 erg s x ) and (c and d) X-ray faint clusters 
(0.6 X 10 44 < L x < 2 X 10 44 erg s" 1 ). X-ray luminosities are from lBohringer et al.l l l2UOOl[200g) . 



reduce the SFA and cannot change galaxy structure while 
the observations require the fraction of bulge-dominated 
galaxies to increase significantly as the clustercentric ra- 
dius decreases. Therefore, even though this mechanism 
can explain the radial variation of the fraction of galaxies 
that are red and inactive, it can not explain the radial 
variation of concentration or central velocity dispersion. 
It is also pointed out that only loosely bound thin clouds 
can be swept out of the galaxy and molecular clouds ar e 
mostly unaffected by this interaction (Quilis c t al.ll20dof ). 
Furthermore, we have shown in Figure [TU] that it is ac- 
tually the galaxy-galaxy interaction that mainly controls 
the color and SFA of the cluster late-type galaxies rather 
than the hot cluster gas. 
6. Viscous stripping / Thermal evaporation 
The cold gas in a late-type galaxy moving through the 
hot intracluster gas can be stripped off by momentum 
transfer betw een the cold disk gas and hot cluster gas 
tiNulsenlll982D. or be evaporat ed by the thermal conduc- 
tion ( Cowie fc Songailal Il977t ). These mechanisms can 



be activated when a galaxy fall into the hot cluster gas 
pool, and can give rise to a characteristic feature in the 
radial variation of galaxy SFA as observed in this work. 
However, they share the same problem the ram-pressure 
stripping mechanism has in that they can not modify 
galaxy structure. 

Furthermore, the hot cluster gas cannot explain the 
environmental dependence of color and SFA of late-type 
galaxies either. If the hot cluster gas predominantly con- 
trols the color and SFA of cluster late types, the con- 
tours in Figure [10] should be vertical independently of 
R n , which is not true. 



7. Strangulation 

According to the current understanding of galaxy evo- 
lution, spiral galaxies maintain SFA by accreting gas 
from their hot gas reservoir. It has been proposed that a 
spiral galaxy fell into the hot intracluster gas can lose 
its hot halo gas and stop supplying gas into its disk 
(jLarson et al.lll980t iBekki et al.ll200ai . Then the galaxy 
just consume its already existing cold gas to form stars, 
and the SFA halts after some time when all cold gas is 
consumed. Compared to the ram pressure stripping, this 
mechanism will quench the SFA in late-type galaxies in a 
delayed action. As being a pure hydrodynamic process, 
it also has the problem that the structural parameters of 
galaxies can not be changed. This mechanism also pre- 
dicts the contours to be vertical (dominance of cluster 
influence) for the early-type fraction, color, and SFA pa- 
rameters in Figures [6] and 1 101 which is not supported by 
observations. 

To inspect galaxy properties for their dependence on 
hot cluster gas in more details, we studied the early-type 
fraction in the R n -R space separately for X-ray bright 
(L x > 2 x 10 44 erg s" 1 ) and faint (0.6 x 10 44 < L x < 
2 x 10 44 erg s _1 ) Abell clusters. Figure HUl shows that 
the magnitude of Je is not particularly higher for X-ray 
bright clusters. This can be another evidence against the 
idea that the hot intracluster gas plays the main role in 
morphology transformation of late-type galaxies in mas- 
sive clusters. 

To summarize, the interaction responsible for the 
MRR/MDR in clusters can be either galaxy-galaxy or 
galaxy-cluster interactions from the point of view of 
the actor, or either gravitational (mass-mass) or hy- 
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drodynamic (gas-gas) interactions in terms of physi- 
cal process. Figure [T2] suggests that the galaxy-cluster 
tidal interactions are responsible for structural and kine- 
matic changes of cluster late-type galaxies toward early 
morphological type. Late types seem to become more 
spheroidal through such interactions. Figure [10] indi- 
cates that the hydrodynamic interactions with early- 
type neighbors are responsible for reddening and SFA 
quenching of cluster late-type galaxies. Therefore, the 
morphology transformation of late-type galaxies seems 
to take place in clusters through galaxy-galaxy hydro- 
dynamic interactions and galaxy-cluster/galaxy-galaxy 
gravitational interactions. However, it is reasonable to 
accept that all above mechanisms are contributing to the 
MRR/MDR in clusters to some extent and a particular 
process cannot fully account for all observational aspects. 

5. CONCLUSIONS 

We have studied the dependence of various galaxy 
properties on the clustcrccntric radius and the environ- 
ment attributed to the nearest neighbor galaxy using the 
SDSS galaxies associated with the Abell galaxy clusters. 
Our major findings are as follows. 

1. There exists a characteristic scale where the galaxy 
properties such as morphology, color gradient, and struc- 
tural parameters suddenly start to depend on the clus- 
tercentric radius at fixed nearest neighbor environment. 

2. The characteristic scale depends on galaxy lumi- 
nosity; the faint galaxies with —17.0 > M r > —19.0 has 
the scale at ~ 3r2oo.ci while the scale is ~ r2oo,d for the 
brighter galaxies with —20.5 > M r > —22.5. 

3. The hydrodynamic interactions with nearby early- 
type galaxies seem to be the main drive to quenching star 
formation activity of late-type galaxies in clusters. We 
do not find evidence that the hot cluster gas is the main 
drive down to the clustercentric radius of ~ 0.1x2oo,ci- 

4. The interaction with the cluster potential and ha- 
rassment are the viable mechanisms that can account 
for the clustercentric radius dependence of the structural 
and internal kinematics parameters. 

Existence of the characteristic scale means that the 
local galaxy number density is not responsible for the 
MDR in cluster because the local density is a smooth 
function of the clustercentric radius and has no disconti- 
nuity in general. The MDR appears working in clusters 
and also in the field because of the statistical correla- 
tion between the local density and the nearest neigh- 
bor distance. What is really working in clusters is the 
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morphology-clustercentric radius-neighbor environment 
relation, where the neighbor environment means both 
neighbor morphology and the mass density attributed to 
the neighbor. 

According to the third and fourth conclusions, the mor- 
phology transformation of the late-type galaxies in clus- 
ters is not due to a single mechanism and has been taken 
place through both hydrodynamic and gravitational pro- 
cesses. 
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APPENDIX 

We calculate the cluster velocity dispersion by using the cluster member galaxies after excluding the interlopers. 
To find the interlopers among the cluster galaxies, we compute 6 for each g alaxy, which indicates the lo cal deviations 
from the systemic velocity (v sys ) and dispersion (a c \ a n) of the entire cluster ([Dressier fc Shectmanl .1988). It is defined 

by 

= 2"" [( w local - V sys ) + (ciocal - CTcl.all) ] I (1) 
°cl,aU 

1 /2 

where N nn is the number of the nearest galaxies that defines the local environment, taken to be iV ga i ' in this study. 
A^gai is the total number of member galaxies in the cluster. The nearest galaxies are those located closest to the galaxy 
on the sky. vi oca i and <7i oca i are systematic velocity and its dispersion estimated from N nn nearest galaxies, respectively. 
We use the galaxies with S < 3.0 to compute the cluster velocity dispersion. 
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